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\ We evaluate the Cronin effect in pA collisions at the CERN LHC and at 

. RHIC in the framework of Glauber-Eikonal model of initial state multipar- 

' ton interactions. Taking carefully into account all kinematical constraints 

■^j- , of each multi-parton interaction process we obtain a softening of the spec- 

f^S ' trum of produced partons, improving in this way the agreement of the 

model with the recent measurements of ir° production in d + Au collisions 
at y/s = 200 AGeV. 

<jj' PACS numbers: 24.85+p, 11.80.La, 25.75.-q 

1. Introduction 

X 

5_j ■ Of particular interest in the context of hadron-nucleus collisions is the 

study of the Cronin effect, which consists of the nuclear modification of 
the transverse momentum spectrum of hadrons with respect to what is 
to be expected from a naive superposition of nucleon-nucleon collisions. 
Because of large scale of exchanged momenta the problem can be approached 
with pQCD methods, while the target complex structure effects may be 
controlled by changing energy and atomic mass number. This effect can 
be explained in terms of multiple scattering of the projectile partons, being 
this latter induced by the high density of the nuclear target (see Ref. [Q 
for a review of theoretical models). With the adoption of the Glauber 
prescription of factorization of the overall many-parton S matrix, which 
can be expressed as a convolution of elementary partonic S matrices, and 
with the decoupling of longitudinal and transverse degrees of freedom, the 
longitudinal component of the incoming projectile parton is conserved and 
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all multiparton interactions of the projectile can consequently be summed 
in the following analytical formula for the inclusive transverse spectrum: 



and Ta(P) is the usual nuclear thickness, as a function of the hadr on- nucleus 
impact parameter b, G(x) the parton number density of the projectile, as 
a function of the fractional momentum x, pt the transverse momentum of 
the final observed parton. The previous quantity can be expressed in terms 
of the dipole-nucleus hard cross section 5^ ar( ^, which originates from the 
square of the scattering amplitude and depends from the transverse size r 
of dipole: 



where do p hard j d 2 pt is the pQCD parton-nucleon cross section, this latter de- 
pending from the parton number density of the nucleus and the elementary 
parton-parton cross section, which includes also the kinematical constraints. 
The infrared divergences deriving from pQCD are regularized with a cut- 
off po; nevertheless, as unitarity is explicitly implemented, the degree of 
infrared singularity of the cross section is reduced from an inverse power 
to a power of a logarithm of the cut-off. It is important to notice that in 
the low pt limit unitarity produces a suppression of the integrated parton 
yield and a random walk of parton to higher pt, recovering in this way the 
local isotropy in transverse space of the black disk limit, which is maximally 
broken in the lowest order impulse approximation. 



As Glauber-Eikonal model does not account for energy conservation, the 
spectrum is therefore shifted towards larger transverse momenta, this effect 
being emphasized as the number of rescattering grows. After implementing 
kinematical constraints exactly, for a given final state, an increased energy 
is needed for the initial projectile partons, and, as the structure functions 
are singular in the small x limit, the initial parton flux is sizably reduced: 
the multi-scattering series of Eq.(^Q) cannot be resummed anymore and the 
expansion in the number of rescatterings is needed: 




(1) 



where 





2. Transverse spectrum expansion 



da drW da^ da^ 



(2) 



d?bdxd?pt d?bdxd?pt d 2 bdxd 2 pt d 2 bdxd 2 pt 



+ ... 
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where the first term of expansion is the single scattering term (the projectile 
parton interacts with a single parton of the target and viceversa) , while the 
other contributions represent the rescattering terms (the projectile interacts 
with i = 2, 3, . . . target partons). Most of the spectrum is nevertheless well 
reproduced by the first three terms of the expansion Ref . |5j . Recalling that 
in the Glauber-Eikonal model the three body process can be expressed as 
a convolution of two on shell two body interactions, it is possible to recon- 
struct the whole kinematics keeping as independent variables the outgoing 
and incoming fractional momenta, the exchanged transverse momenta and 
implementing energy-conservation and mass shell conditions; the following 
j-scattering contribution to the transverse spectrum of i-parton species is 
obtained: 

d 2 bdyd 2 p t ~ Ta ^J II dx 'k ^(?i> • • • > <&) ^iv* x 'v ■ ■ ■ > x j> ?i> • • • Qj) 

x x fi/ p (x,Qfct) fA(x'i,Qfct) • • • fA(x'j, Qfct), 

where A^\qi, . . . , qj) are the subtractive terms deriving from unitarity im- 
plementation and cr 3 (y, x[, . . . , x'j\ q%, . . . qj) are the j-scattering cross sec- 
tions with exact kinematics (see Ref. (31 for a detailed explanation); higher 
order effects in the elementary interactions are accounted by multiplying 
the lowest order expressions in a s by the factor kf ct , while the infrared 
divergences are regularized by a cut-off pq. 



3. Numerical Results 

To study the effect at the LHC we consider the case of production of 
minijets in a forward calorimeter (77 € [2.4,4]) at two different center of 
mass energies in the hadron-nucleon cm. system \fs = 5.5 , 8.8 ATeV . Our 
results are plotted in Fig^ where we compare the spectrum obtained by 
the exact implementation of energy conservation in the multiple interactions 
(solid line), with the approximate kinematics results given by the first three 
terms of expansion of Eq.(|Tj). As an effect of the exact implementation 
of kinematics the spectrum of outgoing particles is shifted toward lower 
transverse momenta; the entity of such a suppression is of the order of 40- 
50% at pt ~ 15 GeV, it is still about 30-38% for higher transverse momenta 
Pt ~ 30 GeV. The triple scattering approximation, used to evaluate the 
spectrum, breaks down at pt < 9 GeV at y/s = 8.8 TeV: by increasing the 
center of mass energy the density of target partons grows rapidly and the 
contribution of higher order rescatterings cannot be neglected any more at 
Pt < 9 GeV (right panel of Fig.l). As the energy is lowered to y/s = 5.5 TeV, 
the region of numerical instability is shifted to the region pt < 1 GeV (left 
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Fig. 1. Transverse momentum spectrum of partons produced in p + Pb collisions at 
^/s = 5.5 , 8.8 ATeV and -q g [2.4, 4], using p = 2 GeV, and kf actor = 2, factoriza- 
tion and renormalization scales Qf c t, Qm equal to the regularized transverse mass 

m t = VhT 
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panel of Fig.l). For a comparison with recent measurements of the Cronin 
effect in dAu — * ttqX at RHIC, we consider the following expression for the 
inclusive ttq spectrum 



da 



frag 



cPqt dy h d?b 



E 



da; 



d 2 pt dy d 2 b 



®Di_> h (Q 2 F ), 



(3) 



where Di^h{QF) are the fragmentation functions at the fragmentation scale 
Qp. At lower energy, = 200 AGeV, we follow Ref. A. in the evaluation 
of the cross section of d + Au — > ir°X, using Qf c t = Q rn = Qf = ^ , 
and the values po = 1-0 GeV and kf act = 1.04. With these choices the 
effects of higher order corrections in a s are minimized and the inclusive 
cross section of ttq production in pp collisions at the same cm. energy is 
reproduced without smearing with the intrinsic hf. the resulting Cronin 
ratio in d + Au — > tt°X is hence a parameter-free prediction of the model. 
By using the leading order K-K-P fragmentation functions at y = and 
b = bdAu = 5.7 fm, which is the estimated average impact parameter of the 
experiment, we evaluate 
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Fig. 2. Left Panel: Cronin ratio in d + Au — > ir° X collisions at y/s — 200 GeV . 
The solid line refers to the case of exact kinematics and the dashed line to the 
case of approximate kinematics at y = 0. The data are form Ref.jSJ- Right Panel: 
Centrality dependence of Cronin ratio with and without energy loss implementation 
at y = 0. Comparison with experimental data presented by PHENIX. 

In the left panel of FigEl we compare our result (continuous line) with the 
experimental data of Phenix collaboration Ref. 5 and with the results using 
approximate kinematics (dashed lines) . Because of the small rapidity values 
of the observed 7To, in this case the corrections induced by exact kinematics 
are important in the region of smaller transverse momenta. The dependence 
of the effect on the impact parameter is shown in the right panel of Fig|2J 
where our calculation (continuous line) is compared with preliminary data 
presented by PHENIX at the DNP fall meeting [5], and with the standard 
Glauber-Eikonal calculation (dashed line). As a consequence of exact imple- 
mentation of kinematical constraints a systematical reduction of the Cronin 
curve is observed, improving in this way the agreement with experimental 
data. 

As it can be seen from the behavior of Cronin ratio in the forward rapidity 
region at r\ = 3.2 FigEl the quenching of the spectrum, due to the energy 
lost by the projectile in the multiple collision of the projectile, is now siz- 
ably increased (notice that the vertical scale is now different respect to the 
previous of Fig. EJ); actually this has to be attributed to the fact that as 
larger become the rapidity values as larger becomes the average number of 
rescatterings. The expectation is nevertheless that the Cronin curve should 
exceed one for pt > 2 GeV; as we can see, from a comparison of our re- 
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Fig. 3. Cronin ratio in d + Au — > tt° X collisions at yfs — 200 GeV . The solid line 
refers to the case of exact kinematics and the dashed line to the case of approximate 
kinematics at r\ — 3.2. The data are form Ref. [7] 

suits with experimental data of BRAHMS Collaboration [7j, for low pt and 
forward values of pseudorapidity, the behavior of the Cronin curve can be 
described with a not so bad approximation also in the context of Eikonal 
dynamics. 

4. Conclusions 

In high energy proton-nucleus collision the Cronin effect is successfully 
described by Glauber-Eikonal model, which allows the implementation of 
unitarity constraints. Working out the leading terms in the expansion in 
multiple parton collisions, we have taken into account all kinematical con- 
straints exactly obtaining an improved agreement with the Cronin ratio 
measured by PHENIX in dAu — ► ttqX at y/s = 200 GeV. Results at larger 
pseudo-rapidities indicate a reduction of the Cronin ratio, which can be 
partially explained by simple Eikonal dynamics. 

REFERENCES 
[1] A. Accardi, [arXiv:hep-ph/0212148 . 

[2] A. Accardi and D. Treleani, Phys. Rev. D 64, 116004 (2001) 
|arXiv:hep-p h/010630"6|. 

[3] E. Cattaruzza and D . Treleani, Phys. Rev. D 69, 094006 (2004) 
|arXiv:hep-ph/0401067| . 



Correct- Proceeding printed on February 2, 2008 



7 



[4] A. Accardi and M Gyulassy, Phys. Lett. B 586, 244-253 (2004) 
|arXiv:nucl-t"h /0308029| . 

[5] S. S. Adler et al. [PHEN IX Collaborationl , Phys. Rev. Lett. 91,072303 (2003) 
|arXiv:nucl-ex /0306021 . 

[6] http:/ /www. phcnix.bnl.gov/conferences. html 

[7] |http://www4.rcf~n l.gov/brahms/WWW/presentations.htmll 

K.Hagel,TAMU : "Forward Physics With BRAHMS at RHIC", CCAST 
2004, Bejing, Aug 10-14, 2004. 



